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Abstract 

We consider the p-adic random walk model in a potential, which can be viewed as a generalization of 
p-adic random walk models used for description of protein conformational dynamics. This model is based on 
the Kolmogorov-Feller equations for the distribution function defined on the field of p-adic numbers in which 
the probability of transitions per unit time depends on ultrametric distance between the transition points 
as well as on function of potential violating the spatial homogeneity of a random process. This equation, 
which will be called the equation of p-adic random walk in a potential, is equivalent to the equation of p-adic 
random walk with modified measure and reaction source. With a special choice of a power-law potential 
the last equation is shown to have an exact analytic solution. We find the analytic solution of the Cauchy 
problem for such equation with an initial condition, whose support lies in the ring of integer p-adic numbers. 
We also examine the asymptotic behaviour of the distribution function for large times. It is shown that in 
the limit t —>■ oo the distribution function tends to the equilibrium solution according to the law, which is 
bounded from above and below by power laws with the same exponent. Our principal conclusion is that the 
introduction of a potential in the ultrametric model of conformational dynamics of protein conserves the 
power-law behaviour of relaxation curves for large times. 


1 Introduction 

An ultrametric space is a metric space M with metric d{xi,X 2 ), xi, X 2 C M satisfying the strong triangle 
inequality 

Va:i, X2X3 e M : d{xi,X2) < max{(i(a;i, 0 : 3 ), d{x2,X3)}. (1) 

A metric d{xi,X 2 ) satisfying ([T]) is called an ultrametric. A classical example of an ultrametric space is the field 
of p-adic numbers Qp [Dill [3]. The p-adic numbers were introduced in 1889 by K. Hesel and were proved widely 
useful in algebraic geometry, theory of numbers and representation theory. 

For a long time the ultrametric approach have been found useful in solving various problems in the field 
of classification of objects and information processing of data arrays [3]. In the last 30 years the ultrametric 
analysis had had a recent reemergence—it received a great impetus from the pioneering works of researches 
from the scientific school of Academician V. S. Vladimirov, whose efforts were later taken up by a number of 
researches from different scientific schools (for an overview, consult, for example, [S]). This relative new research 
field is now known as the “p-adic and ultrametric mathematical physics” and is blessed by a number of books 
and an immense number of research articles in the area of p-adic analysis, p-adic mathematical physics and 
their applications to modeling in various areas of physics, biology, computer science, sociology, physiology, etc. 
(see, for example, [Sl[5] and the references given therein). 

First ultrametric models in physics emerged first in the 1970s in the theory of spin systems with disorder 
niziiais]. It was found that if a system has a large number of multiscale “internal contradictions” (frustrations), 
then it may reach an equilibrium in hierarchically nested regions of the phase space, the nesting level increases 
with decreasing temperature. Here, the ratio of scales of phase domains satisfies the strong triangle inequality, 
and thus the low-temperature spin states are found to be ultrametrically correlated. Almost immediately 
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after the emergence of the ultranietric model of spin glass, a conjecture was made about the ultrametricity of 
the space of conformational states of a protein molecule [iniiii]. It should be noted, however, that a literal 
transfer of the ultrametric description of phase states of the spin glass and conformational states of protein 
molecule was triggered by a then-popular analogy between the low-temperature proteins and glasses [12] . As 
in the case of glasses, the atomic mobility of polymer globules is accompanied, in particular, by multiscale 
constraints generating highly intersected energy landscapes. In this sense, proteins resemble glasses, some of 
their thermodynamical properties are the same (for example, the low-temperature behaviour of heat capacities). 
But, as distinct from disordered globular structures, proteins have a unique special feature—they are capable of 
performing precise operations with objects of atomic scale, the physical reasons behind it remaining unclear. 

A systematic treatment of the conformational dynamics of proteins based on an ultrametric model is given 
in a series of papers by the authors [EllIllllSlIlelllilllHlIiniEolEl]. In these papers, the conformational 
dynamics of protein was described as follows. The conformational state of a protein was considered as a quasi¬ 
equilibrium macrostate, which includes many microstates of a macromolecule. A conformational state which is 
understood in exactly this sense is also called a “basin”. Further, there is a great number of bonds on mechanical 
degrees of freedom of atoms in a macromolecule. As a result, there are many constraints on the flexibility of 
a macromolecule at various scales of the configuration space of degrees of freedom of the entire macromolecule. 
It is assumed that a random walk of a protein over the configuration space is effected as follows. The most 
probable are such motions of atoms in a macromolecule which take place in small regions of the configuration 
space (small-scale rearrangements of atoms). In order that there would be a rearrangement of atoms on a large 
scale involving a much larger region of the configuration space, it is necessary that rearrangements of smaller 
scales would result in a configuration of a macromolecule from which a passage to a different region of the 
configuration space would be enabled by the available bonds. The probability of such an event is much smaller 
than the transition probability between arbitrary small-scale configurations. Hence, the characteristic times of 
probability transitions between configurations are much smaller on more fine scales that those of probability 
transitions between configurations on more coarse scales of the configuration space of a macromolecule. In 
accordance with this approach, one may conventionally single out an increasing family of scales (characteristic 
sizes of the regions) of the conhguration space of a macromolecule. Correspondingly, the characteristic times 
of a macromolecule to stay in regions of the configuration space of various scales increase with increasing scale. 
From such a model of a random walk of a protein over regions of the configuration space (conformations or 
basins) of various scales one may catch the hierarchy of transitions and the hierarchy of characteristic times of 
the transitions. The indexation of such a hierarchy enables one to introduce in a natural way an ultrametric 
distance between conformational states. This model features a number of simplifying assumptions. Namely, 
it is assumed that the distance graph is a Kelley tree and that all conformations are homogeneous in energy. 
These assumptions are fairly restrictive, but nevertheless, they enable one to formalize such a model by a p- 
adic random process. In doing so it is assumed that the probabilities of transitions per unit time through the 
activation barriers (that is, transitions between conformational marostates) obey the Arrhenius relation. Such 
a model, which looks simple at first sight, has enabled one to adequately describe two principal experiments: 
the experiment on spectral diffusion in proteins m and the experiment on binding kinetics of (70 by myoglobin 
[Min]. We note that in these two experiments, the fluctuation dynamics of a protein is described in a unified 
way. 

The further development of this model consists in its generalization by cancelling one or several simplifying 
assumptions. One of the ways towards it is to reject the assumption on the homogeneity of conformations 
relative to energies. This can be achieved by introducing a potential in the model. 

The statement of the problem of ultrametric diffusion in a potential well is discussed in the literature on p-adic 
mathematical physics starting from the 2000s. Here, there is a principal difficulty regarding the introduction 
of a potential in the equation of an ultrametric random walk. It is known (see, for example, |22| 1 that the 
right-hand side of the equation for the distribution function of an arbitrary time-homogeneous Markov process 
contains three terms. The first term describes the drift of a trajectory and is related with the potential force, 
the second term describes the diffusion processes, and the third term describes the jump processes. In the 
case when the trajectories of a random process are everywhere continuous, the term responsible for the jump 
processes vanishes and the equation becomes the Fokker-Planck equation, which contains two terms—the term 
responsible for the potential force acting on a walking particle and the diffusion term. Since a p-adic random 
walk is supported by functions that are constant everywhere except at discontinuities of the first kind, it follows 
that the first two terms, including the term containing the potential, vanish. Containing only the third term. 
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the equation now becomes the equation, which is conventionally called the Kolmogorov-Feller equation or the 
master equation. In this case, it is natural to introduce a potential in the form of some additional structure on 
an ultrametric space; this additional structure should be incorporated in the probability of transition per unit 
time, which in turn enters the master equation. 

One way to introduce potential in the model of a p-adic random walk is to consider the master equation for 
the probability density of the form 

df{x,t) [u{y)f{y,t)-u{x)f{x,t) 

“ J -‘ ’ 

Qp 


where Qp is field of p-adic numbers, dpy is the Haar measure on Qp, a is positive parameter interpreted as the 
inverse temperature, u (x) is positive definite function on Qp (potential function). We assume that the function 
u (x) belongs to the class Wp, for some j3 < a i.e. satisfies the properties: 1) |u (x) | < C (l + Ixl^); 2) there 
exists n, such that u{x + x') = u (x) for all \x'\p < p"". The Cauchy problem for equation ([5]) is determined 
in the class of functions f{x,t), which are continuously differentiable in t, and as a function of x G Qp belong 
to the class class Wj for 7 < a — /3 uniformly with respect to t. In this work, equation @ will be called the 
equation of p-adic random walk in a potential u (x). 

At present, methods of analytical solutions of equation @ in potential of arbitrary form u (x) are unknown. 
Nevertheless, it turns out that for a definite form of a power-law potential this equation has a precise analytic 
solution. 

This article is devoted to the analytical solution of the Cauchy problem for equation ([5]) with the initial 
condition on the ring of integer p-adic numbers Zp for the specific choice of the potential function u (x) depending 
only on p-adic norm |x|p of x € Qp as the power law of a special form. This model can be viewed as a 
generalization of the models used in previous works for description of protein conformational dynamics mui 
[13 [la [HI [n [m uni m • In Section 2 we show that, for such special potentials, the equation of p-adic random 
walk in a potential can be reduced to the equation of p-adic random walk with modified measure and reaction 
sources, whose support lies in Zp. The last equation will be the subject of research of the last two sections. In 
Section 3 we examine the equation of a p-adic random walk with arbitrary measure and a source in Zp. Using 
the basis of eigenfunctions of the Vladimirov operator with modified measure, as was constructed in our recent 
paper [23] , we find an analytic solution of the Cauchy problem for such an equation with an initial condition 
in Zp. In Section 4 we particularize the solution of this equation in the case when the measure is induced 
by a potential of special form, which was considered in Section 2. In Section 5 we examine the asymptotic 
behaviour of the so-obtained solution for the distribution function for large times. Namely, we establish that 
the distribution function tends, as t ^ 00 , to the stationary solution according to a law that is bounded from 
above and below by power laws with the same exponent. Our main conclusion is that the introduction of 
a potential in the ultrametric model of conformational dynamics of protein conserves the power-law behaviour 
of the relaxation curves for large times. In Appendixes A, B, C and D we state and prove a number of assertions 
that are used for solving the principal equation of the model and in obtaining an asymptotic estimate of the 
solution thus obtained for large times. 


2 Particular case of a power-law potential 

We consider the equation of p-adic random walk in the potential ([2]) when the function u (x) is in the class Wp 
and has the form 


where 


u (x) = aU (|x|p) -I- b\x\^ (1 - U (|x|p)), 0 < /3 < o. 


U (|x|p) 


1 , |x|p < 1 , 
0 , |x|p > 1 . 


(3) 


We denote (Qp,dpx) C Wq the class of locally constant functions of functions square-integrable with 
respect to p-adic Haar measure dpX. Also we denote Lf ^ {Qp,u (x) dpx) the class of locally constant functions 
of functions square-integrable with respect to measure u (x) dpX. Let us state the Cauchy problem for ([2]) with 
the potential ([3]) and the initial condition 
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( 4 ) 


/ {x, 0) = /o (x) 

in the class (Qp,dpx). The following result holds. 

Lemma 1. Let constants a, b, a and jd satisfy the conditions 


fd = a — 1. 

Then the solution of the Cauchy problem for equation 


df {x, t) 
dt 



{y) dpy + xn (|a;|p) f{x, t), 

\x - y\p^ 


(5) 

( 6 ) 

(7) 


with the initial condition m (here A = bVp (a) Tp (—a) and rp(—a) = — - - - — is p-adic gamma function) in 

the class L(^ (Qp,u{x)dpx) is the solution of the Cauchy problem for equation fUj in the class Lf ^ {Q.p,dpx) 
with the initial condition 0- 

Proof. First of all, note that Lf ^ (QpjW (a:) dpx) C Lf ^ {Qp,dpx). Next, we show that under conditions ([5|) and 
([^ the equation ([2|) can be transformed to the form Q. The equation ([2|) can be rewritten in the form of the 
equation of p-adic random walk with the modified measure u (y) dpy and reaction source as follows 


df (a:, t) 
dt 



(y) dpy + J (x) f(x, t). 

\x y\p 


( 8 ) 


where 


J (x) = Pp (—a) {x). 

is the pseudodifferential Vladimirov operator [5] 


(9) 


D°^lp{x) = 


1 


_ f ‘p{y)-p jx) 

Ppi-a) J \x-y\p+^ 


dpy. 


To calculate D°‘u (cc) we find the p-adic Fourier image u (k) of the function u (x). Using the relations 

/ dpxLl (|x|p) X (kx) = ft (|fc|p) 


and 


we find that 


J dpx\x\^ (1 - U (|a;|p)) x {kx) = 

Qp 

= Tp (a/3 + 1) \k\;d-^ - (I^Ip) - Tp M + 1) (1 - i\k\p)) 


= aTl (|fc|p) -I- b 


u (k) 


J dpxu (x) X (kx) = 


Qp 


1 — p ^ 

1 — p“^“^ 


U(|fc|p)+Fp ild + l)\k\;d-i 
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We have 


D°‘u{x)= j dpk\k\pu{k)x{—kx). 


and hence 


D<^u {x) = (^a - D^n (|a;|p) + bTp (a) (|x|p) 

In view of conditions ([S]) and (ED, we find that 

(x) = bFp (a) (|a;|p). 

As a result, 

j(*)= 


\x - y\ 


a+i dy = Fp (-a) D°‘u (x) = bFp (a) Fp (-a) Ft (|a;|p) = XFl (|a;|p) 


and so (ED assumes the form ([7]) , which completes the proof of Lemma [T] 


□ 


Remark 1. Note that in the framework of assumptions (ED and (ED of Lemma[T]we have A > 0 for a > 1, A < 0 
for a < 1, and A = 0 for a = 1. Moreover, for a < 1 we have /3 < 0, hence we will study only the case a > 1. 


Remark 2. For a > 2 the equation has the normalized to unity stationary solution 


(a:) = 


1 - 1 - p- 


1 -p 


1—p“i p““2 _ xy yi—p 

Later we will show that the equilibrium solution of equation (ED 

[x) = lim / (a;, t) 


— fl(|a;|p) + |a;|p“+ (1 - (|a;|p)) 


j{eq) — i U*) I a >2, 
Do, 1 < a 


< 2 . 


( 10 ) 


1 «> 2 ^ 

1 0, 1 < a 


< 2 . 


Solution of the Cauchy problem for the equation (ED in the class (Qp, u (x) dpx) can be found analytically. 
Below we shall find and explore an analytical solution of this problem for the homogeneous initial condition in 


Zp. 


3 The general solution of the equation of p-adic random walk with 
general modified measure u (x) dpX and reaction source in 


In this section we solve the Cauchy problem of equation 


df (a), t) 
dt 



{y) dpy + Xn (|a:|p) f{x, t) 

\x-y\p 


( 11 ) 


with the initial condition 

f {x,0) = n{\x\p) ( 12 ) 

and an arbitrary function u (x). In finding the solution we shall use a number of auxiliary results, which will 
be formulated and proved in Appendix A. 
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The solution of the Cauchy problem (llll) -(fT ^ will be searched in the subclass of locally constant functions 
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of the space Lnorm (Qp! u (x) dpx) (see Appendix A). Passing to the Laplace images, we obtain 

sf{x,s) = n{\x\p) + D^f^^)f{x,s) + Xf2{\x\p)f{x,s). 

Employing Theorem [T] from Appendix A, we expand the function f (x, s) in basis (I39|) . 


(13) 


/ (a;, s)= ^ /i (s) (j)i (x). 


(14) 


An application of Theorems [2] and |3] and Corollary [T] from Appendix A shows that 


sfk (S) = VoV;_\ (^1 - ' - (l h is) + 


+A VoV-_\ ( 1 - 


14-1 


-1 


1 - 3-5^ + 2 ^''“^ 


Vk 






i—k 


K;- 




V,, 


fk (s) + 


-2XVoK-Ml-^jfk(s) + 


+ ^VoV-_j 1- 


Vk-l\ 2 
Vk 


J2fds)V-J (l- 


V-l \ 2 


for fc > 1 and 


OO OO 


\i—k 


i—k-\-l i—k-\-l 


^fk (s) = - (l -p-(“+i)) i j (s) + A ^ ^ fk is) (16) 


for fc < 1. 

Equation (|15l) can be rewritten as 


(15) 


fk (s) = gk is) + Xgk (s) ^ f^ is) E,_2 f 1 - 


P.-l \2 




(17) 


where 


Setting 


gk is) = VoV,_j 1 - 


Vk-] 

Vk 


y s+(i-p-(“+i)) - 


-AEol 4 -\(l- 


Vk-] 

Vk 


-1 


1 - 3-5^ + 2 


Vk 


i—k 


-X^V,V-\ (l - + 2 AEo14--i (1 - 




Vk-1 

Vk 


Fis)=Y.his)V~_l (l- 


1 ^: 


2-1 


(18) 


Gis) = ^gds)V,_l ( 1 - 


2=1 


1 ^. 


we write equation (1171) in the form 


fk is) = gk is) (1 + XF is)). 


(19) 
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-i / \ 2 

Multiplying (ITO)) by V^_i I 1-) and summing in fc, we find that 




and hence, 


Substituting (1^ into (fTOl) . this gives 


i^(s) = G(s) + AG(s)F(s), 

l-AG(s)' 

7 _ 9k (s) 

1 - AG (s) ■ 

From equation (1161) we find that 

fk (s) = 0, fc < 1. 

In view of (ED and (|22[) . function (|14p assumes the form 

/(».») = 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


Expression (I23p is the general solution of the equation of p-adic random walk with source in Zp with the initial 
distribution in Laplace images. 

4 Solution in the case of a modified measure induced by a power-law 
potential 

In this section we particularize the solution (1231) with the following function u (x): 


(l*lp) + l^lp Ml - {\x\p))^ ■ 

Calculating the volumes of balls (HOI) with respect to the measure u (x) dpX, this gives 


y = 6(i-p-i)^' 


,a(z+l) 


p“ — 1 

Here, the basis functions dMl) assume the form 

_ 1 

(j)i (x) = (^6(1- p-i) ^ (f2 (|a;|pp"*+^) - p““G (|a;|pp"*)) . 

Calculating functions (fTOl) . we obtain 

— ak(2 


a (.) = (i,p“ (1 - p--))’ 


The function G (s) reads as 




(p“ - l)p- 


^ s — bp~'^p°^Tp{—a)' 
Substituting (1271) and (1051) into (IQTl) . we find that 


fkis) = {bp^^ {l-p ^))2 


P 


—akl2 


1 


s — bp~'^p°'Tp{—a) 




(p“ - l)p- 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


s — bp-^p°^Tp{—a) 
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Next, in view of we have 




p “'"p" (f? (|a;|pp ^+^)-p “J7(|x|pp '=)) 


fc=i (s - 6p-'=p“rp(-a)) 1 - AX;“i 


(p“ - l)p- 


s — hp~^p°‘Tp{—a) 


Here, the Laplace image of the probability of finding the trajectory of a random process in the support of the 
initial distribution is as follows: 


/ (x, s) dpX = ^ 


P 


-ak /„a 


(p“ - 1 ) 


=1 (s - hp-^p°‘Tp{-a)) 1 - A X;i=i 


(p“ - l)p- 


s — bp ^p°‘Tp{—a) 


Gis) 


l-AG(s) 


= Fis). 


(30) 


5 Asymptotic behaviour of the distribution function as t oo 

We first examine the asymptotic behaviour as t —>■ oo of the probability of finding the trajectory in the support 
of the initial distribution, which we denote by S (t) and which is the Laplace preimage of function (1301) . To 
this aim we need to find the poles and residues of the function G (s) (see (1^51) 1. Clearly, the function G (s) is 
positive at zero, and besides, 

G ( 0 ) =_= i 

6p“rp(-a)^^ A 


The function G (s) vanishes at infinity 
The derivative of G (s) is negative 


lim G (s) = 0. 


ds 


g(») = -E 


(p“ - l)p- 


^ (s - 6p-*p“rp(-a)) 


< 0 , 


and hence the function G (s) is decreasing. The function G (s) has simple poles at the points Sk = bp~^p'^Tp{—a) < 
0, which correspond to the zeros of the function F (s). To the left from the poles Sk = bp~^p°‘Tp{—a) < 0 the 
function G (s) is negative, and to the right, is positive. The points of intersection of the function G (s) with 

the line G = — correspond to the poles of the function F (s). Besides, the function F (s) has an evident pole 
A ^ 

at the point s = 0, at which the function G (s) assumes the value —. Clearly, all the poles of the function F (s) 

A 

are negative, except for the pole s = 0. The function G (s) vanishes at the inhnity, and hence so is the function 
Fis): 

lim F (s) = 0. 

s—^oo 

Let s = —Ai, i € Z_|_, be the poles of the function F (s). Note that the function F(s) is not meromorphic 
(the point s = 0 is a point of condensation of simple poles). Nevertheless, one may show (see [TH] for details) 
that the function F(s) can be expanded in simple poles. 




bk 


k=l 


Afe 


Clearly, Ai satisfy the condition 

-6p-(*+i)p“rp(-a) < A, < -bp-^p^Vp{-a). 


We write Ai as 


A. = Vrp(-a) (1 + (5.) 










where Si is a small positive error satisfying the condition 


From the equation for the poles, 


it follows that 


0 < < p — 1. 

1 


( 31 ) 


G(-Afc) = 


A’ 


ftri — k — 1 


P~ 


p-k-1 _\_ p-k-1^^ _ p-i pa-1 _ \ 


(32) 


Equation (15^ enables one to obtain an asymptotic estimate for Sk as k ^ oo (see Appendix B, formula (1551) 1. 
The residues at the poles s = —Xi of the function F (s) are as follows: 


-2 




p“ - 1 

- 1 )" 


E- - - 

^ (p-fc-i _p-^) 


(33) 


Equation (j33p in view of asymptotic estimate for Sk enables one to obtain an asymptotic estimate for bk as 
A: —)• oo (see Appendix C, formula ([55)) b This estimate can be written as 


(a) (1 + 5^,2 (fc-2 _ 1)) (1 + o (1)) < 

<bk< G,nax (a) (1 + 5„,2 { k -^ - 1)) (1+0 (1)) , 

where o ( 1 ) is an infinitesimal quantity as fc —>■ oo and 

' (p“ - 1) (1 (1 -p-“+2)^ 


(34) 


c„,„(o) = I (jFzI) 

(P- 1 ) 


p2-a (pa-1 _ ly 


2 ’ 


Gmax (cr) — (1 + O (1)) < 


- 1) si,^^ 

(p2-a (pa-1 _ 

(p“ - 1) (1 -p-“+2)^ 
p2-a (pa-1 _ ly 

(l -p-^) (l -p-^)' 

ip - 1 )^ 

(p° - 1 ) 

^ p2-a (p“-l — 1) 


a > 2, 
a = 2, 

1 < a < 2, 

a > 2, 


2 ’ 


, a = 2, 

1 < a < 2. 


The function S (t) is as follows: 

OO 

S{t) = bo-\-^bk exp (-Xkt). 

k^l 

Using estimates for bj- and for Xk one can estimate S{t) as 

bo + Gmin (a) (l + Sa ,2 {k~^ - l)) cxp {-wp~'"p~^t) (1 + 0 ( 1 )) < 

k 

< S{t) < 6 o + Gmax(a)E^’”'“~^'^ (1 + < 5 a ,2 “ l)) exp {-wp~'^p~^t) (1 + o(l)), 


(35) 


where 


w = — 6 p“rp(—a). 
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An application of Lemma 2] from Appendix D gives, as f > c», 


^0 + A (a) t ^1^1 + Sa ,2 ^+ ^ 

< S{t) < bo + B (a) t 1“ + Sa ,2 ^ ^ ~ ) (1 + 0 (1)) 


( 36 ) 


where A (a) and B (a) are defined by 


A (a) = Cinin (a) P ^'(Inp) ^ {wp “ 'r(|a-2|), 
B (a) = Cmax (a)p“'““^'(Inp)”^ (wp“^) ^'r(|a-2|) 


and o (1) is an infinitesimal quantity as t —> c». 

We note that 

(p“ - 1 ) (p““^ - l) 

ho= ' 


(p“-i - 1 )^ 


0 , 


a > 2, 

a < 2, 


is nothing else but the equilibrium solution at the point a; = 0, i.e. 

lim / (0, t) = lim S (t) = bo, 

t—^oo t—>-oc 

which for a > 2 concides with the normalized stationary solution /+*) (a;) at the point x = 0 

lim S (t) = bo- 

t—^CO 

Let us discuss the asymptotic behaviour of the distribution function / (a;, t) as t —>■ oo in the domain \x\p > 1. 
For the Laplace image of / {x, t) we have 


f{x,s) = 


K (s, \x\p) + G (s) 
l-AG(s) 


(37) 


where the function K (s, |a;|p) is defined as 


cn ™ —a 


^ (s, |a;|p) = -- 


p mp 


In |x|p/ Inp 


(s-6|a;|p V“rp(-a)) 


- E 


i=l 


p-ai (pa _ 

(s - 5p-*p“rp(-a))’ 


and the function G (s) is defined by formula (1281) . Clearly, the function (l37l) has a pole at the point s = 0, at 
which the residue bo (x) is the equilibrium solution of equation ([7]), which for a > 2 concides with the normalized 
stationary solution /+*) (a:). From (1571) it follows that the function f {x,s) also has poles Afe, which coincide 
with the poles of the function F (s). The residues bk (x) at these poles are as follows: 


bk (kip) = lim (s + Afc) 


—^ — 


K {s,\x\p) + G {s) 
(l-AG(s)) • 


(38) 


Formula 

where 


can be rewritten as 


bk (kip) = (p“ -1) 


bk (kip) — i,bk (kip) + ^) bk. 


p“klp“ 


(p_fc_l p-fe-14 _ \x\p 


In |a:|p/ Inp 

E 


p-az (pa _ 


i=l 


(p-/c-l _|_p-fe-l^^ -P“®) 


and bk is the residue at the pole s = —Xk of the function F (s). For any fixed kipi the function hk (kip) has 
finite nonzero limit as k ^ 00 . Hence, as k ^ 00 the behaviour of the residues bk (kip) bk are equal up 
to a factor depending on kip- For this reason, if t ^ 00 and for any x G Qp the function / {x,t) tends to the 
equilibrium solution in accordance with the same law as the function S (t). 
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6 Conclusion 


In the present paper we consider the equation of p-adic random walk in a potential, which has the form ([5]) . It 
is shown that in the case when the function of potential has a power form ([3]) and in the particular case when 
the parameters are related by ® and m, the equation ([5]) can be solved analytically. We found the solution 
to the Cauchy problem of this equation with the initial condition in Zp and examined its asymptotic behavior 
at large times. Moreover, we prove that, as f ^ oo, the distribution function tends to the stationary solution 
according to a law bounded from above and below by power laws with the same exponent. The analysis of 
numerical solutions of the p-adic random walk in a potential |24| . which we not present here, shows that, for 
other values of the parameter /3 near the value /3 = a — 1, the power-law behaviour of the distribution function 
is conserved for large times. 

The results obtained in this work provide a basis for the belief that the introduction of a power-law poten¬ 
tial in the p-adic model of the conformational dynamics of protein conserves the power-law behaviour of the 
relaxation curves for large times. This fact may serve as a justification of the conclusion that in a model that 
is capable of describing experiments in spectral diffusion in proteins and experiments in binding kinetics of CO 
by myoglobin (in which the relaxation has power-law behaviour), the introduction of a power-law potential may 
not be necessary. We admit, however, that in more precise experiments on the conformational dynamics of 
protein the potential may somehow exhibit itself. It is quite possible that the potential may play a certain role 
in models of prototypes of molecular nanomachines that are based on processes of p-adic random walk. One 
of such models, which was considered in m, was based on a process of p-adic random walk without potential, 
but on the compact set C Qp. Such a model can be physically interpreted as a model on Qp in a potential 
u(x) of a special form, which is constant for x G and is infinite for x ^ B^. 
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Appendix A. Some auxiliary results 

Let u[x) is positive dpX-locally integrable function on Qp. We denote by L^Q,,^(Qp, u{x)dpx) the class of 
functions on Qp of the form / (x) = f (|a;|p) with finite integral 

J dpyu (y) \f{y)f . 

Qp 

Theorem 1, The family of functions 

_ 1 

(x) = (^U_i (^1 -(^l7(|a;|pp-*+^) --5^l7(|x|pp-*)^ , (39) 

where f G Z and 

Vi = f m (y) dpyQ (|p|pP“*) , (40) 

Qp 

forms an orthonormal basis for (Qpi (®) dpx). 

Proof. That the system of functions (1331) is orthonormal 

J 4>i (x) 4>j (x) m (x) dpX = Sij 

Qp 


11 




is verified directly. Let us show that (|5^ is a basis for ^ i'^) dpX). Clearly, the countable family of 

the supports of all p-adic spheres 

= {x e Qp : |x|p = p’'} 


of the form 


^ {\x\p-p1 


1, X S Sr, 
0, X ^ Sr 


forms a basis for L'^orm (Qpj i^) dpx). We have d (|x|p — p’’) = 17 (|x|pp“’') — 17 (|x|pp“’'+^), and so it suffices 
to show that, for any r, the function 17 (|x|pp“’') can be expanded in functions (j39|l . We claim that the norm of 
the following function is zero. Indeed, 


C (Ixipp-’') - W £ V~J ( 1 - 


i—r-i-1 


V, 


4>i (x) 


j 17 (|x|pp ^)nn{x)dpX— 

Qp 

i 

-2Vr ^ ( 1 - y 17 (|x|pp“’') (/)i (x) m (x) dpX+ 

< 0 p 

{(j)i (x))^ TO (x) dpX = 


i—r+1 


+ v ^ E ^-1 1 - 


i—r+l 


V,-i 


= W - 2W E ^-1 ( 1 - 


i=r+l 






+ Vr^ ( 1 - 


i—r-i-2 


V.JVrll- 




V,-i\2 


V, 


+ 




= Vr-Vr^ £ V-l - £ = 0 




i—r+2 


As a result, we have 


c (Ixipp-’-) = w E ^-1 ( 1 - 

completing the proof of the theorem. 




C-l \2 




(x) 


(41) 

□ 


Any function / (x) from L’^^rm (Qpj (x) dpx) can be expanded as 


OO 

fix)= > 

i— — co 


where 


/i = y h (x) f (x) TO (x) dpX. 


We note that functions (15^ have the property 


y (j)i (x) TO (x) 


dpX = 0 . 


The following result holds [23]. 
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Theorem 2. Functions (I39|l are eigenfunctions of the operator that generalizes the Vladimirov operator in 
Ll (Qp, rn (x) dpx) 


D’^{x)f{x) = J dpym{y) 


f{y)-f jx) 


(42) 


with the eigenvalues 

OO 

= - (l 

j^i 


(43) 


Lemma 2. The following formula holds 


/ 


(fi (x) (fj (x) 4>k (x) m (x) dpX = 


where 


— dij5jkVj._-^ I 1 


Vk-i 

Vk 




Vk 


Vk^ 




Vk- 


Vk 


+ SikSi<:jyj_i ( 1 — 






+ 


+SjkSj^iV^_^ ( 1 


V ^-1 \ 2 




dj<i — 



j < i, 
j > *• 


The lemma is is proved by direct calculation. 

Theorem 3. Let f{x),g (x) G (Qp: u (x) dpx), 

OO OO 

f{x)= Mi (x), g(x) = g^(j)i (x). 

2 = —OO i — — oo 


Then 


where 


OO 

f{x)g{x)= ^ hkcfkix), 

k— — oo 


hk = fkgkVkJy 1 - 


X4-1 

Vk 


1 - 3-^^ + 


Vk 


Vk 



Vk- 


Vk 


1 


fidi “t” fk 


i—k-\-l 


1 



-f- 


OO I 

+gk 

i—k-\-l 



1 

2 
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Proof. We write 


OO 


OO OO 


/ (x) g{x)= Y {x) cj)j (x) = Y ^k(t>k {x), 


where 


2 = —OO J — — 00 


OO OO p 

hk= Y X! 


k— — oo 


Z = —OO J — — 00 

'S.p 

Now the theorem follows from Lemma [5J □ 

The following result is a direct corollary to Theorem [3] taking into account the expansion formula 

^ Vo 


^{\x\p) = Yy.^^ 


i=0 


Corollary 1. Let f (x) € (Qpi ^ (x) dpx) and f (x) = J2T=-oo fkf’k (x). Then the production of functions 

fi (|x|p) / (x) is expanded in the basis (j39l) as 

OO 

n{\x\p)f{x)= Y Ck4>k{x), 


where 


Ck = fkVoV^-_\ ( 1 - 


Vfc-] 

14 


k— — c 


-1 


1 - 3-5^ +2^''“^ 


Vk 




+h E y«yr-\ (i - -y^) - (i - ■^) 


+ y^vti 


for k > 1 and 


for k < 1. 


Vk 


Ck = fk Y ^0^-1 (1 - 

i—k-\-l 


V 


V-1 

V 


(44) 


(45) 


Appendix B. Estimate of Sk for large k 

To estimate Sk for large k it is convenient to write equation (|32ll as 


p-{a-l)kp-a+l ^ p-{a-l)i 

Sk ^ 


^ l_p-fc-l+*(l+4) S l_p-i+fc+l(l + 4) pa-1- 

2 — 1 2 — 

In what follows we shall require the following result. 

Lemma 3. // 0 < x < 1, then 


p—aipk+l 


(46) 


1 


If Q < X < y <1, then 


1 — X 

1 

1 — X 


> 1 + X. 


< 1 + cx, 


(47) 

(48) 


where c = 


1 - 2 / 
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Proof. 


1 




1 — X 


1 — X 


1 — X 


1 — X 1 — X 


□ 


Consider the case a > 2. By (SZl) and ((151) and Lemma [5] with a 7 ^ 2 we estimate the terms on the right 
of (1151) as follows: 

_ 1 )^ ^ 

i: 1 (1+4) > ( 1 +(1+ V) = 

2 = 1 ^ ^ 2=1 


E 


p 


1 _ p-(“-i)fc 

pO-i _ X 
— (a — 1)2 


+ (1 + Sk) 


1 _ p-(«- 2 )fe 

pa-2 _ I 


-P 


^l_p-»-.«(l + 4) 


< (1+ 


l_p-{a-l)k I ^ I _ p-ia-2)k 

pa-1 _ I p _ X p““2 _ I P ’ 


Moreover, the following estimate holds: 


°° ^-ai k+l °° „—ai k+l 

E P P _ _ P P 

1 _ „-i+fc+i n^ x-B-i(i 


p—{a — l)kp—2a+l 


i=/c+2 ^ + + (l-P “)(l + 4)' 


- E —^ 

2 —^ 1 — n— 


p—aipk+1 


°° „-(a-l)fc„-2a+l 

< _ p-a^ k+l ^ _p _ p _ 

X _p-*+fc+i (x + 4) ^ p p X 

i=k+2 ^ i=k+2 


■P 


Applying these estimates, we have 


p-{a-l)kp-a+l ^ p-{a-l){k+l) (X + 4) (l - ^-(“"2)'=) 


4 


>- 


p-{oL-l)k p-2a+l 


1 — p““+l 


pa -2 _ X 


(1 -p 1) (1 -p “) (1+4)’ 


(49) 


-(a-l)fe -a+1 -(a-l)fe 1 -1-A, 1 _ r,-(a-2)fe -(a-l)/c -2a+l 

- -+-< r + , L, , P - ^ , + —■ (50) 

Ofe p“^ — 1 p— lp“^ — 1 1 —p“ 

Next we need the following assumption on the behaviour of (5^ as fc —>■ oo. We make the ansatz that Sk —> 0 
as k —>■ oo and solve these inequalities with this ansatz. Then, in the asymptotic limit k oo, the last two 
inequalities may be written as 

(1 - p-i) (1 - p-“+2) p-(“-2)fe (X + o (1)) < 4 < (1 - p-“+') p-(“-2)'= (1 + O (1)), a > 2. (51) 

Consider the case a = 2. In this case we estimate the terms on the right of (ll 6 )) as 


E 


p 




X-p-'=-i+4i + 4) , . 

2=1 2=1 


>y]p-*(i+p-'=-i+4i+4)) = 


-j _^ 

= 1 + (1 + '^fc) 

pct 1 _ 1 


-k-1 


V_— 

Z_^ 1 — n~^~ 


— (a— 1)2 


<E+(l + Ayir+"'+''- 


^ l_p-i-i+i(l + 4) - i_p 

1 -+ +i±4t.p-‘, 


p““^ — 1 p — 1 
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p-aipk+l 


p-2zpfc+l 


p >^p 3 


_ >- V _ - _ =_ 


-Sr 


p—aipk-\-l 


< 


. — 7 ?“*+^+^ (1 + Sk) 

i=fc+2 ^ ^ ^ i=fc+2 


- ^ p ‘"p'"' ^ = 


^—2i^k-\-l 


p kp ^ 
1 — 


Using these estimates, we may write 


—k —1 _U _1 

P ^ , >-^^ + (l + 4)A:p-'=-i 


„-k 


4 


1 -p- 


(p- 1) (p2 - 1) (1 + 4) 


— i’ —1 

p Up i. 


< -- 


, 1 + 4, _fe p V ^ 

■ H- —kp - 


4 " p“-i-l p-1 l-p-2- 

We again make make the ansatz that —>■ 0 as fc —^ oo. With this ansatz the last two inequalities read as 

(l — p~^) k~^ (1 + o ( 1 )) <4 < fc~Ml + o (1)) I a > 2 . (52) 

Consider the case 1 < a < 2. One may show that inequalities (|4^ and ([50l) with the ansatz Jfc ^ 0 as 
fc —> oo lead to a contradiction. Hence, in this case we make a different ansatz that Sk —>■ const as k ^ oo. This 
gives us the inequalities 


P 


,-(a-l)fcp-a+l ^_(a_l)(fc+l) (1 + 4 ) _ p-(a- 2 )fc^ 


> - 


4 1 ~ p p“ ^ — 1 


p—(a — l)kp—2a+l 

(1 -p-1) (1 -p"“) (1 + 4)’ 


{a-l)kp-a+l p-(a-l)k 1 + 4 1 — p-(“-2)/c ^ ^ p-{a-l)kp-2a+l 


Sk 

which may be written as 


pal_l'l_pl pOi 2 _ 


1 -p- 


Sl+{1- 


p — 1 p — 1 

1 _p-a+l “ (1 -p- 1 ) - 1 ) 


) 4 - (p-“+" - 1 ) < 0 , 


<5^+11- 


(1 - 4 - (1 - p-‘) - 1 ) 

4 - (p-“+' - 1 ) = 0 


- 1 ) > 0 . 


We denote by Jmax the maximal positive root of the equation 

p-a +2 _ ^ — 1 


V 1 - P"“+^ (1 - P”^) (P“ - 1) 

and denote by i5min the maximal positive root of the equation 


4 - (i - ^ 


- 1 ) = 0 . 


Then 


^min (1 + 0 (1)) < 4 < (^max (1 + 0 (1)) . 

One may show that limQ ,^2 <5max = 0, lima^i (5max = oo and limQ ,^2 ^min = 0, lima^i i5niin = oo. 

Combining these results, we write down the lower and upper estimates for Sk as k ^ oo for different a as 

dmin ^ Sk ^ D max 5 (53) 


where 


'(l-p-^) (l-p-“+ 2 )p-(“- 2 )^ a > 2 , 
U>mi„ = (l + o(l))<j (l-p-i)fc-i, a = 2 , 

(5mini 1 + Cr + 2, 


D„ 


■(l_p-a+ 2 )p-(a- 2 )fe, q, > 3 , 

= (l + o(l))U-i, a = 2 , 


1 < a < 2. 
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Appendix C. Estimate of bk for large k 

Let us estimate the residues as fc —>■ oo. To this aim we write (1331) as 


bk = 


- 1 

(p“-i - 1 )" 


V- 


+p-^-^6k -p~") 


E 


(p-fc-l ^p-k-l§^ -p“*) 


p 2 -ap-(a- 2 )fe 


Taking into account that in the case a > 2 we have i5fe ^ 0 as fc —> oo, we write 


(54) 


bk = 


p°‘ - 1 

(P“-' - 1 )' 


E—^ 


E 


(p-fc-l _ p-i-j 


p 2 -ap-(a- 2 )fe 


n -1 


( 1 + 0 ( 1 )). 


Let us estimate bk from the above. First, we consider the case a > 2. Using (1551) . we find that 

-1 


bk < 


p°‘ - 1 

(P“-' - 1 )' 


E- — -2+ E 


p-ai p2-ap{a-2)k 


.fcl {P ^ ^-P , i^k+2 

Summing and neglecting the terms that vanish for fc —> oo, this gives 

fr)“ - 1 m 

bk < 


{p->^-^f (l-p-“+2)^ 


( 1 + 0 ( 1 )). 


p 2 a ^pOL 1 _ 2^^2 


For a = 2 we similarly have 

In the case 1 < a < 2 we have 
- 1 


(p- 1 ) 


bk < 


(p“-i - 1 )' 


p-a _ p-a(k+l) 


P 


-a(k+ 2 ) 


p‘^-ap-{a-2)k 


which gives 


(1 — p~°‘) — p~^) (1 — p~°‘) {p~^~^) 

0 ( 1 )). 


e; 


-1 


( 1 + 0 ( 1 )), 


bk < p-(2-^)fc (1 + , 

p2-a (pa-1 _ 1)2^ 


Let us now estimate bk from below. For a > 2, we have 
- 1 


bk > 

(p“-i - 1 )^ 
which gives 


P 

-k-l 


.S {p-'^-'^ - p-^Y 


E 

i—k -\-2 


p 


p2-ap{a-2)k 


(p-fe-l_p-i)2 [1 — p 1)^(1—p a+2(2 


(1 + 0 ( 1 )). 


( p “ - 1 ) (1 - p -')^ (1 - ,1 + „, 1 )) ^ 


p2-a (p“-l — 1) 


In the case a = 2, we have 


bk > 


E - 1 ) ,.-2 


(P - 1 )' 


For 1 < a < 2, we get 

p“ - 1 


bk > 


(pa-i _ 1)2 


p-“l _ p-a(fc+l) 


A )-2 (1 + 0 ( 1 )). 


p-akp-2ap-ai 


{1 — P °‘) {p ^ ^ — p (1 — p “) (p ^ 1 — p *^ ^) 


p 2 - 0 !p-(a- 2 )fc 

^max 


(1 + 0 ( 1 )), 
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and hence. 


b > _(p! - l)'^max -(2-o)fc ^ ^ _ 


Collecting the above inequality, we write the lower and upper estimates for 6^ as fc —>■ oo for different a as 
follows: 

a {a, fc) (1 + o (1)) <bk <b {a, fc) (1 + o (1)), (55) 


where 


(p“ —l)(l—p p “+2)^ 


: (a, fc) = (1 + o (1)) < 


p2 “ (p“ ^ — 1) 

(P' - 1 ) /,-2 

(P-If ’ 

(p“-l)(5max (2-aU 

I p2-“ (p“-l - 1)2 


a = 2, 

1 < a < 2, 


b {a, fc) = (1 + o (1)) < 


(p“ —l)(l—p “+2)2 ^ 


p2-a (pa-1 — 1) 

(l-p-2) (1-p-l)^ 

IP-if 

(P°‘ - 1) ^min (2-c)k 
[p2-“(p“-l-l)2^ 


p-(a-2)k^ a > 2 , 


a = 2, 

1 < a <2. 


Appendix D. Asymptotic estimate of the series 

Here, we shall obtain the estimate for the series 

OO 

Si{t) = ^ exp (^—b~^t) , t>0, a>l, 6>1, fc>0 (56) 

i=l 


for large t. 

Lemma 4. For series (1561) we have the following asymptotic estimate as t —^ oo : 

a~^ (In6)^”^ t“ln& (Int)”^ F ((1 + o (1)) < Si{t) < a (In6)^”^ t“TnT (Int)”^ F ((1 + o (1)), 

\ In 6 / V In 0 / 

where F (z) is the Gamma function and o(l) is an infinitesimal quantity as t —> 00 . 

Proof. Instead of series (1551) . we consider the series 


<S' 2 (t) = ^* *exp(—6 ’'t) = Si{t) — a ^exp(—6 ^t) . 


We note that —ra “ is a decreasing function and exp (—6 ^t) is an increasing function of x. Hence, on the 
x'^ 

interval i < x < i + 1 we have 


-^a ^ exp (^—b < a * exp (—6 < 


(x- 1)' 


-a 


-(x-l) 


exp {-b-^t) . 


(57) 


Integrating (1571) in x from i to i + 1, we find that 
ri+i 


^ ( —a exp (—b dx < a *exp(—6 < a / - - TTk^ “exp (—6 ^t) dx. (58) 

Ji ^ ^ ^ Ji 
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Next, summing (1581) in i from 2 to oo, this gives 


where 


(t) < 52 (t) < 5niax(i); 

OO 

5min(0 = a~^ J exp dx, 

^ exp {-b-^t) dx. 


5max(^) — ^ J' ' 


(x — 1)^ 

2 

In ((59l) and (|60)l we change the variable of integration as = y and b~^t = y, respectively. We have 


(59) 

(60) 


b-H 


Smin{t) = a ^ (ln6) Inf) 


Inc 


-k 


Int — Inw \ inn_, 

-+ yinb ^exp{-y)dy, 


5max(0 = a(ln&) t Inb 

If fc > 0, then, as t oo, 


b-^t 

_i _ina f f\nt — \ny 


Inb 


-k 


Inc 


-1 y^nb exp{-y)dy. 


Sniin{t)=a ^(ln6)^ inb (Int) '"r ( ) (l + o(l)) 


Inc 


-k - 


In a 
In b 


5max(t) = a (Info)^ Inf) (In t) '‘r ( ) (1 + o(l)), 


In a 


\-k - 


In a 


which gives the assertion of the lemma. 


□ 
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